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Abstract 

Improving previous calculations, we compute the J/f) tt —► charmed mesons cross section using 
QCD sum rules. Our sum rules for the J/if it —> D D*, D D*, D* D* and D D hadronic matrix 
elements are constructed by using vaccum-pion correlation functions, and we work up to twist-4 in 
the soft-pion limit. Our results suggest that, using meson exchange models is perfectly acceptable, 
provided that they include form factors and that they respect chiral symmetry. After doing a 
thermal average we get (a nJ ^v) ~ 0.3 mb at T = 150 MeV. 
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I. INTRODUCTION 


For a long time charmonium suppression has been considered as one of the best signatures 
of quark gluon plasma (QGP) formation [lj. Recently this belief was questioned by some 
works. Detailed simulations [Q] of a population of c — c pairs traversing the plasma suggested 
that, given the large number of such pairs, the recombination effect of the pairs into char- 
monium Coulomb bound states is non-negligible and can even lead to an enhancement of 
J/'ijj production. This conclusion received support from the calculations of [|[]. Taking the 
existing calculations seriously, it is no longer clear that an overall suppression of the number 
of J/-0’s will be a signature of QGP. A more complex pattern can emerge, with suppression 
in some regions of the phase space and enhancement in others ||, |5j. Whatever the new 
QGP signature (involving charm) turns out to be, it is necessary to understand better the 
J/^j dissociation mechanism by collisions with comoving hadrons. 

Since there is no direct experimental information on J/-0 absorption cross sections by 
hadrons, several theoretical approaches have been proposed to estimate their values. In 
order to elaborate a theoretical description of the phenomenon, we have first to choose the 
relevant degrees of freedom. Already at this point no consensus has been found. Some 
approaches were based on charm quark-antiquark dipoles interacting with the gluons of 
a larger (hadron target) dipole |?|. ||] or quark exchange between two (hadronic) bags 
0,0. whereas other works used the meson exchange mechanism [11], [12| 0QO^- i" 


this case it is not easy to decide in favor of quarks or hadrons because we are dealing with 
charm quark bound states, which are small and massive enough to make perturbation theory 
meaningful, but not small enough to make non-perturbative effects negligible. Charmonium 
is “the borderguard of the mysterious border of perturbative world of quarks and gluons 
and the non-perturbative world of hadrons” [[17] . 


In principle, different approaches apply to different energy regimes and we might think 
that at lower energies we can use quark-inter change models S0 or meson exchange models 
Tl|, [T2|, [l3|, n~4|. [I5i [TTill and, at higher energies we can use perturbative QCD || [7. §. 


However, even at low energies, the short distance aspects may become dominant and spoil a 
non-perturbative description. In a similar way, non-perturbative effects may be important 
even at very high energies JT8[ . 

Inspite of the difficulties, some progress has been achieved. This can be best realized 
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if we compare our knowledge on the subject today with what we knew a few years ago, 
described by B. Mueller (in 1999) [[19] as “...the state of the theory of interactions between 
J/ip and light hadrons is embarrassing. Only three serious calculations exist (after more 
than 10 years of intense discussion about this issue!) and their results differ by at least two 
orders of magnitude in the relevant energy range. There is a lot to do for those who would 
like to make a serious contribution to an important topic”. In the subsequent three years 
about 30 papers on this subject appeared and now the situation is much better, at least in 
what concerns the determination of the order of magnitude, which, as it will be discussed 
below, in the case of the J/ip pion interaction, is determined to be 1 < < 10 mb in 

the energy region close to the open charm production threshold. 

One of the main things that we have learned is the near-threshold behavior of crj/^_ 
This is quite relevant because in a hadron gas at temperatures of 100 — 300 MeV most of 
the J/'ij) — Ti interactions occur at relatively low energies, barely sufficient to dissociate the 
charmonium state. In some calculations a rapid growth of the cross sections with the energy 
was found ||T2| , p~3] [L4|] . This behavior was criticized and considered to be incompatible with 
empirical information extracted from J photoproduction [j20|. This criticism, however, 
made use of the vector meson dominance hypothesis (VDM), which, in the case of charm, 
is rather questionable [J2TJ. The introduction of form factors in the effective Lagrangian 
approach, while reducing the order of magnitude of the cross section, did not change this 
fast growing trend around the threshold. Later, again in the context of meson exchange 
models, it was established |L6] that the correct implementation of chiral symmetry prevents 
the cross section from rising steeply around the threshold. In a different approach, with 
QCD sum rules (QCDSR) 


J, the behavior found in |16] was confirmed and in the present 
work, with improved QCDSR, we confirm again the smooth threshold behavior. We, thus, 
believe that this question has been answered. 

Another, phenomenologically less important, but conceptually interesting issue is the en¬ 
ergy dependence in the region far from threshold. Results obtained with the non-relativistic 
quark model |9] indicated a rapidly falling cross section. This behavior is due to the gaussian 
tail of the quark wave functions used in the quark exchange model. This result of the quark 
model approach could be mimicked within chiral meson Lagrangian approaches with the 
introduction of y/s dependent form factors [[T^, |23| . For J/ijj — N interactions, it was found 
in ref. H] that this behaviour depends ultimately on the gluon distribution in the proton at 
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low x. In the case of J/ijj — N, for certain parametrizations of the gluon density one could 
find a falling trend for the cross section but no definite conclusion could yet be drawn. 

If the J/ijj is treated as an ordinary hadron, its cross section for interaction with any 
other ordinary hadron must increase smoothly at higher energies, in much the same way as 
the proton-proton or pion-proton cross sections. The underlying reason is the increasing role 
played by perturbative QCD dynamics and the manifestation of the partonic nature of all 
hadrons. Among the existing calculations, no one is strictly valid at y/s ~ 20 GeV, except 
the one of ref. |R|, which is designed to work at very high energies and which gives, for the 
.//^-nucleon cross section the value crj /^-n — 5 mb. This number can be considered as a 
guide for J/-0 — tt cross section in the high energy regime. It should, however, be pointed 
out, that the calculation of ref. ]T8| is based on a purely nonperturbative QCD approach. 
The inclusion of a perturbative contribution will add to the quoted value and will have a 
larger weight at higher energies. A similar conclusion was reached in |^4[. In the traditional 
short distance QCD approach the cross section grows monotonically [6|, [7. ||. 

As a side-product, the theoretical effort to estimate the charmonium-hadron cross section 

27], within the framework of QCD sum rules, of 


motivated a series of calculations [25, 


form factors and coupling constants involving charmed hadrons, that may be relevant also 
to other problems in hadron physics. 

In this work we improve the calculation done in ref. |22] by considering sum rules based 
on a three-point function with a pion. We work up to twist-4, which allows us to study the 
convergence of the OPE expansion. Since the method of the QCDSR uses QCD explicitly, 
we believe that our work brings a significant progress to this important topic. 

The paper is organized as follows: in the next section we review the method of QCD 
sum rules, giving special emphasis to the QCD side. In section III we present some formulas 
for the computation of open charm production amplitudes and in section IV we give our 
numerical results. Finally some concluding remarks are given in section V. 


II. THE METHOD 


In the QCDSR approach [[2^, [29], the short range perturbative QCD is extended by an 


OPE expansion of the correlator, giving a series in inverse powers of the squared momentum 
with Wilson coefficients. The convergence at low momentum is improved by using a Borel 
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transform. The coefficients involve universal quark and gluon condensates. The quark-based 
calculation of a given correlator is equated to the same correlator, calculated using hadronic 
degress of freedom via a dispersion relation, giving sum rules from which a hadronic quantity 
can be estimated. 

Let us start with the general vaccum-pion correlation function: 


n Al34 — 


d 4 x d 4 y e ip2 ' y e ip3 ' x {Q\T{j 3 (x)j^)j^(y)}\'K(p 1 )) , (1) 


with the currents given by fjf = cy^c, 73 = wr 3 c and j 4 = cT 4 <i. p \, p 2 , Pz and p 4 are 
the four-momenta of the mesons 7 r, J/^, M 3 and M 4 respectively, and T 3 and T 4 denote 
specific gamma matrices corresponding to the process envolving the mesons M 3 and M 4 . 
For instance, for the process J/ijj n —> D D* we will have ^ = 7 ^ and T 4 = 775 . The 
advantage of this approach as compared with the 4-point calculation in ref. [Q, is that we 
can consider more terms in the OPE expansion of the correlation function in Eq. ( 0 ) and, 
therefore, check the “convergence” of the OPE expansion. 

Following ref. |30[, we can rewrite Eq. ([[]) as: 
d A k 


n M34 — — 


(2tt 


Tr[S ac (p 3 - k)^ lx S cb (p 3 ^p 2 - fc)r 4 T> afe (fc,pi)r 3 ] , 


where 


Sabip) = i 


. j> + m c 


p 2 — m 2 


Jab 


( 2 ) 


(3) 


is the free c-quark propagator, and D ab (k,p) denotes the quark-antiquark component with 
a pion, which can be separated into three pieces depending on the Dirac matrices involved 
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D ab (k,p) = S ab 775-4 + 7a75^“ + 7 5 o- q/ 3 C'" /j ] , (4) 

with a, b and c being color indices. The three invariant functions of k, p are defined by 


A(k,p) = y 2J d ^ X e * k ' X ( 0 \ d ( x ) ir Y5 u (°)W( p)) , 

B a (k,p ) = ^Jd 4 x e^(0|d»7 Q 75«(0)|7r(p)) , 
C af3 (k,p) = -^f d A x e ik - x {0\d{x)a a ^ 5 u{0)\Tr(p)) . 


(5) 


Using the soft-pion theorem, PCAC and working at the order 0(p^p u ) we get up to 
twist-4 


A (k,p) = 


( 2 tt ) 4 (qq) 

12 U 


-2 + ip ai — 


d 


idk, 


Oil 


1 d d 

+ 3 PaiPa2 idk ai idk, 


OL2 


<5 (4> (fc) , 
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B a (k,p) = 


( 2 vr) 
12 


-fn 


iPo 


1 di5 2 f \ 

“1“ 2 PaPai ~gj7 ^ "gg" \^Pol9ol\ol 2 ~ ^P^Qaai J T 


d_ _ d_ 

idk ai idk a2 


ri n . „\ ( 27r ) 4 (99>,„ „ N 

Cap\k,p) \Pa9/3a 1 P(39aa 1 ) 


. a 




i<9fc ai 2 i(9fc ai i<9/c Q2 


<5 ,4) (fc) , 


<5 (4) (fc) , 
( 6 ) 


24 3/, 

where S 2 is defined by the matrix element {^\dg s Q°‘^^pu\'K(j))) = id 2 f n p a , where Q a p = 
e aP orQ aT /2 and Q a0 = t A G a p. 

The additional contributions to the OPE comes from the diagrams where one gluon, 
emitted from the c-quark propagator, is combined with the quark-antiquark component. 
Specifically, the c-quark propagator with one gluon being attached is given by [p9[ 

9 sQo.fi 


[kola ~ kp'Va + + m c )i<7 al 3 } , 


(7) 


2 (k 2 — ml) 2 

Taking the gluon stress tensor into the quark-antiquark component, one can write down the 
correlation function into the form 

d A k 


n^34 — —4 


(2vr)' 


-Tr 


Safi{P3 - k)^f^S(p Z -p 2 -k) 


+S(p 3 - k)'y fl S a /3(p3 Pi fc) )r 4 L> Q/ 3 (fc,pi)r 3 


where we have already contracted the color indices, and where we have defined 

1 


Sopik) = 


[k a 7/3 - kpria + Qt + m c )ia a p\ , 


and 


with 


2 (k 2 — m 2 ) 2 

D a/3 (k,p) = ^a pX E pXa(S (k,p) +Ye a0es F T9 s(k,p) , 

E pXa P(k,p) = ~ Jd 4 x e ik - x (0\d(x) l5 a pX g s G^u\7r(p )), 
F T es(k,p) = ^ J d A x e lk - x ($\d(x)-f T g s Ge 5 u\'K{;p)). 


( 8 ) 


(9) 


( 10 ) 


( 11 ) 


Up to twist-4 and at order 0(p p p u ), the two functions appearing above are given by 
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E pXap = -^f3n (p a p p g xis - P 0 P p g Xa - P a P x g pP + P 0 P x g pa ) (2n)W 4 \k) 


32 

if ) 2 f 

F-m = - -p s g re )(2^‘\k) , 


( 12 ) 


where f 3n is defined by the vacuum-pion matrix element (Q\dg s a a p^§Q al 3 u\Ti{p)) 


6 





















The phenomenological side of the correlation function, 11 ^ 34 , is obtained by the consider¬ 
ation of J/^, M 3 and M 4 state contribution to the matrix element in Eq. ([[]). The hadronic 
amplitudes are defined by the matrix element: 

iM = <^(p 2 ,aOI M 3 (-p 3 ,v) M 4 (-p 4 ,p) 7r(pi)) 

= i Mtfn(jai,p2,P3,P4) 4.IT.IT , (13) 

where /*“ = e*“ for the D* meson and /*“ = 1 for the D meson. 

The phenomenological side of the sum rule can be written as (for the part of the hadronic 
amplitude that will contribute to the cross section) | 22 | : 

j-fphen __ 3 A 4 Mfj ,34 _ 1 /-i^N 

^ 34 (pi - m\) (pi - mi)(pi - ml) 

where h. r. means higher resonances, and where A, is related with the corresponding meson 
decay constant: (T>|j D |0) = -A D = m 2 D f D /m c and (0|j a |T>*) = X D *e a = m D *f D *e a . 


III. HADRONIC AMPLITUDES FOR J/ty vr -*• OPEN CHARM 

The hadronic amplitudes can be written in terms of many different structures. In terms 
of the structures that will contribute to the cross section we can write 

• for the process J/-0 tt —> D D *: 

Mpv = Af D *Pi M Pm + X-2 D *PipP2u + A 3 D *p lu p 3/1 + a + A 3 D *P2uP3p , (15) 


for the process J/'i/j 7r —> D D: 


^-DD ^fia( 3 crPi P3P4 1 


(16) 


for the process J/^ 7 r —> D* D*: 


Mfwp = Af " D * H^p + Jp Vp + A%* D * g^pCpap-yPiP^pl + Af * D * e upa/3P3pPl P 3 

A e vpaPP 3 pPlP2 + Aq* D *€ p Va i3p 3 pP^P2 + A^* D * e^apPlpPiP^ 

A 3 * D * eupapPluPlPA + Ag* D * 6p iU paPi + Af 0 * D *e^ a/3 pi M p 4 p^ + A®l D * epvpaP 2 
Ai2 D 4,,ua0PlpPlP3 + Af 3 ‘ D ‘e^p 3 pPipf ■ (17) 
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with H pvp = (evafagw - epa^g^p^pj + e lipa pp 2v piP^ and J pup = {e vpafj p Ul + e ppoif5 p lv + 
Cfiva 0Plp)P2P3 + tp.ua/3P2pPlPz- 

In Eqs. (H), (|1|) and ([l7|), A t are the parameters that we will evaluate from the sum 
rules. In principle all the independent structures appearing in H pup and J pup would have 
independent parameters A,. However, since in our approach we get exactly the same sum 
rules for all of them, we decided to group them with the same parameters. 

Inserting the results in Eqs. (||) and © into Eqs. (^) and (||) we can write a sum 
rule for each of the structures appearing in Eqs. flT5|), ( |T6j ) and ([17]). To improve the 
matching between the phenomenological and theoretical sides we follow the usual proce¬ 
dure and make a single Borel transformation to all the external momenta taken to be equal: 
— P 2 = —p\ = —p\ = P 2 —> M 2 . The problem of doing a single Borel transformation is the 
fact that terms associated with the pole-continuum transitions are not suppressed fh|. In 
the present case we have two kinds of these transitions: double pole-continuum and single 
pole-continuum. In the limit of similar meson masses it is easy to show that the Borel 
behaviour of the three-pole, double pole-continuum and single pole-continuum contribu¬ 
tions are e~ m M / m 2 /M 4 , e ~ m ^ //M2 /M 2 and e~ - m M / M2 respectively. Therefore, we can single 
out the three-pole contribution from the others by introducing two parameters in the phe¬ 
nomenological side of the sum rule, which will account for the double pole-continuum and 
single pole-continuum contributions. The expressions for all 19 sum rules are given in the 
Appendices A, B and C. 


IV. RESULTS AND DISCUSSION 


The parameter values used in all calculations are m c = 1.37 GeV, m n = 140 MeV, rripj = 
1.87 GeV, m D * = 2.01 GeV, = 3.097 GeV, U = 131-5 MeV, (qq) = -(0.23) 3 GeV 3 , 
ml = 0.8 GeV 2 , <5 2 = 0.2 GeV 2 , f 3w = 0.0035 GeV 2 |k|. For the charmed mesons decay 
constants we use the values from |32| for /d and fp* and the experimental value for /^: 


= 270 MeV, f D = 170 MeV, f D * = 240 MeV. 


(18) 


In ref. |34|] we have analyzed the sum rule for the process J/^ tt —> D D. Here we choose 
to show the sum rule for A^* D * in Eq. (£T|), as an example of the sum rules for the process 
J/V> 7T —> D* D*. 
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FIG. 1: Sum rule for A®* D related to the process J/i/j tt —>■ D* D* as a function of the Borel 
mass. The dots, squares and diamonds give the twist-2, 3 and 4 contributions to the sum rule. 
The triangles give the result from Eq. 0. The solid line give the fit to the QCDSR results. 


In Fig. 1 we show the QCD sum rule results for A±* D * + Af* D * M 2 +B®* D * M A as a function 
of M 2 . The dots, squares and diamonds give the twist-2, 3 and 4 contributions respectively. 
The triangles give the final QCDSR results. We see that the twist-3 and 4 contributions are 
small as compared with the twist-2 contribution, following the same behaviour as the sum 
rule for the process J/rj) tt —> D D given in [34 . In general all the other sum rules are similar 
and contain twist-2, twist-3 and twist-4 contributions corresponding to the first, second, and 
third terms inside the brackets in the right hand side of Eq. m- Only the sum rules for 
up to Af^ D *, A® D * and A® D * do not get the leading twist contribution, and will be 
neglected in the evaluation of the cross section. It is also interesting to notice that if we 


consider only the leading twist contributions we recover the sum rules obtained in ref. |22 
The triangles in Fig. 1 follow almost a straight line in the Borel region 6 < M 2 < 16 GeV 2 . 
This show that the single pole-continuum transitions contribution is small. The value of the 
amplitude Af* 12 * is obtained by the extrapolation of the fit to M 2 = 0 |26| , |27|, |33j . Fitting 
the QCD sum rule results to a quadratic form we get 


A°* d * ~ 10.5 GeV 3 . (19) 

Since we worked in the soft pion limit, A®* D *, as well as all other A, is just a number. All 
particle momenta dependence of the amplitudes is contained in the Dirac structure. 
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In obtaining the results shown in Fig. 1 we have used the numerical values for the meson 
decay constants given in Eq. (0. However, it is also possible to use the respective sum 


rules, as done in [E3J. The two-point sum rules for the meson decay constants are given in 


the appendix D. The behaviour of the results for the hadronic amplitudes does not change 
significantly if we use the two-point sum rules for the meson decay constants instead of the 
numerical values, leading only to a small change in the value of the amplitudes. In Fig. 2 
we show, for a comparison, both results in the case of Af >D . 



FIG. 2: Sum rule for Af * D * related to the process J/i[) it —» D* D* as a function of the Borel mass. 


The dots and squares give the results from Eq. (Cl) when using respectively numerical values or 
the two-point sum rules for the meson decay constants. The solid and dot-dashed lines give the 
fits to the QCDSR results. 


Using the respective sum rules for the meson decay constants we get 


Af* D ' ~ 13.9 GeV 


-3 


( 20 ) 


We will use these two procedures to estimate the errors in our calculation. It is important 
to mention that our results agree completely with the value obtained in |22[ . 


The results for all other sum rules show a similar behaviour and the amplitudes can be 
extracted by the extrapolation of the fit to M 2 = 0. The QCDSR results, evaluated using 
the numerical values for the meson decay constants, as well as the quadratic fits for the 
amplitudes associated with the process J/-0 tt —> D D* are shown in Fig. 3. 
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FIG. 3: Sum rules for Af D * related to the process J/ij) ir —> D D* as a function of the Borel mass. 
The dots, squares, diamonds, triangles up and triangles down give the results from the QCDSR for 
ki D up to AjP D respectively. The QCDSR results for Ai D and A± D can not be distinguished 
at the scale of the figure. The solid, dotted, dashed, long-dashed and dot-dashed lines give the fits 
to the QCDSR results. 

The values for the parameters associated with the process J/ip n —> D D* are given in 
Table I. 


TABLE I: The best fitted values for the parameters associated with the process «//?/> n —> D D*. 


A 

A° D * 

AD D * 

Af B * 

A|p D * 

14 ± 2 GeV -2 

-7.2 ± 0.9 GeV -2 

-58 ±8 GeV" 2 

14.6 ±2.2 

-15.6 ± 2.2 GeV -2 


For the process J/-0 n —> D D we have only one parameter which is given by [[II 


A dd = 13.2 ±1.8 GeV" 3 , 


( 21 ) 


and the 13 parameters associated with the process J/ip n ^ D* D* are given in Table II. 
The errors in all parameters were estimated by the evaluation of the sum rules using the 
numerical values and the two-point QCDSR for the meson decay constants. 

Having the QCD sum rule results for the amplitudes of the three processes J/i^ 7r —> 
D D*, D D, D* D*, given in Eqs. (|T5P , (pb|) and ( |T7|) we can evaluate the cross sections. 
After including isospin factors, the differential cross section for the J/-0 — tt dissociation is 
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TABLE II: The best fitted values for the parameters associated with the process J/ijj n —> D* D*. 


A D*D* 

A°* d * 

a d* D * 

a d* d * 

A mm 

12.2 ± 1.7 GeV 3 

-12.8 ± 1.8 GeV“ 3 

12.5 ±1.7 GeV 3 

-24.6 ±3.4 GeV 3 

9.8 ± 1.6 GeV 3 

a d*d* 

a d*d* 

Af D ‘ 

Ag* D * 

\D*D* 

A 10 

9.7 ± 1.6 GeV" 3 

-13.0 ± 1.8 GeV 3 

-13.8 ± 1.8 GeV 3 

-5.4 ±0.9 GeV 1 

2.5 ±0.2 GeV 3 

Afi^GeV 1 ) 

l\D*D* 

IV 12 

A D*D* 
a 13 



(—5.5 ± 0.5)10 -3 

-0.022 ± 0.002 GeV 

3 0.53 ±0.03 GeV 3 



given by 

da 

1 Vl\4l 2 


(22) 


dt 

967T5D? 

yu/l ^Pz,cm spin 



where p i)Cm is the three-momentum of p\ (or p 2 ) in the center of mass frame (with pi (p 2 ) 
being the four-momentum of the 7r 


2 A(s,m2,mJ) 

Pi,cm = - 


4s 

with A(x, i/, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz, s — (pi + p 2 ) 2 , t = (pi — p 3 ) 2 . 
In Eq. d22j), the sum over the spins of the amplitude squared is given by 


EV = av-H’v hr'- 


spin 


a u! 

P2P2 

ml 


' »v> _ PsP 3 _ 

y 9 

V m n* 


for J/ip Ti —> D D*, with p 3 (p 4 ) being the four-momentum of (£)). 


E I-UP = 


spin 


for J/ip 7T —> £) and 


EV = hr' 


spin 


LL LL 

P2P2 

ml 


LL LL 

P2P2 

ml 


P 3 P 3 


P 4 P 4 


m 


D* 


m 


D* 


(23) 


(24) 


(25) 


(26) 


for J/ip n —> D* D' 

The structures multiplying A{ ,rr and A^ D ' in Eq. (|l^), and in Eq. HD break 

chiral symmetry |16|. To evaluate the effect of breaking chiral symmetry in the process 
J/-0 7T —> D D* + D* D we show, in Fig. 4, the cross section calculated using all structures 
in Eq. (|T5| ) (dashed line) and neglecting A® D * and Aj? D * (solid line). 
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FIG. 4: J/tf} 7r —> D D* + D* D cross section. The solid and dashed lines give the results which 
respect and break chiral symmetry respectively. 

From Fig. 4 we see that the cross section obtained with the amplitude that breaks chiral 
symmetry grows very fast near the threshold. Since this is the energy region where this kind 
of process is probable more likely to happen, it is very important to use models that respect 
chiral symmetry when evaluating the J/'ip — it cross section. As mentioned before, the sum 



FIG. 5: J /ip — n dissociation cross sections for the processes J /ip ir —> D D* + D D* (solid line), 
J /ip it —> D D (dashed line) and J/ip tt —> D* D* (dot-dashed line). 

rules for A ® D *, A^ D * and A^ D * up to A^ D *, do not get the leading order contribution and 
will be neglected when evaluating the cross sections. It is important to keep in mind that, 
since our sum rules were derived in the limit p\ —> 0, we can not extend our results to large 
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values of y/s. For this reason we will limit our calculation to y/s < 4.5 GeV. 

In Fig. 5 we show separately the contributions for each one of the process. Our first con¬ 
clusion is that our results show that, for values of yfs far from the J/-0 tt —> D* D* threshold, 
aj/fn->D*D* — a J/ 4 >n^DD*+DD* > i 11 agreement with the model calculations pre¬ 


sented in [0 but in disagreement with the results obtained with the nonrelativistic quark 
model of [[|, which show that the state D*D has a larger production cross section than 
D*D*. Furthermore, our curves indicate that the cross section grows monotonically with 
the c.m.s. energy but not as fast, near the thresholds, as it does in the calculations in 
Refs. Hl2| , [I3| , pj p| , pj . Again, this behavior is in opposition to ||, where a peak just after 
the threshold followed by continuous decrease in the cross section was found. 



FIG. 6: J/tj) — tt dissociation cross sections from meson exchange model (k| (dot-dashed line), 
quark exchange model (dashed line), short distance QCD §1 (dotted line) and QCD sum 
rules (solid lines). The shaded area give an evaluation of the uncertainties in our calculation. 

In Fig. 6 we show, for comparison, our result for the total cross section for the J/tf} tt 
dissociation (solid lines) and the results from meson exchange model [T2|] obtained with 
a cut-off A = 1 GeV (dot-dashed line), quark exchange model || (dashed line) and short 
distance QCD SI (dotted line). The shaded area in our results give an evaluation of the 
uncertainties in our calculation obtained with the two procedures described above. It is very 
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interesting to notice that bellow the DD* threshold, our result and the results from meson 
exchange and quark exchange models are in a very good agreement. However, as soon as 
the DD* channel is open the cross section obtained with the meson exchange and quark 
exchange models show a very fast grown, as a function of -y/s, as compared with our result. 
As discussed above, this is due to the fact that chiral symmetry is broken in these two model 
calculations. 

The momentum distribution of thermal pions in a hadron gas depends on the effective 
temperature T with an approximate Bose-Einsten distribution. Therefore, in a hadron gas, 
pions collide with the J/iJj at different energies, and the relevant quantity is not the value of 
the cross section at a given energy, but the thermal average of the cross section. The thermal 
average of the cross section is defined by the product of the dissociation cross section and the 
relative velocity of initial state particles, averaged over the energies of the pions: (a nJ ^v), 
and is given by [|T| 




— l£0 


f™dz[z 2 - («i + a 2 ) 2 ][z 2 - («i - a 2 ) 2 ]K 1 (z)a nJ ^(s = z 2 T 2 ) 


(27) 


4a 2 K 2 (ai)a 2 K 2 (a 2 ) 

where cq = rrii/T (i — 1 to 4), £ 0 = max(aq + a 2 , a 3 + aq) and K t is the modified Bessel 


function. 



FIG. 7: Thermal average of «//?/> dissociation cross section by pions as a function of temperature 
T. The shaded area give an evaluation of the uncertainties in our calculation. 

As shown in Fig. 7, (a nJ ^v) increases with the temperature. Since the J/if) dissociation 
by a pion requires energetic pions to overcome the energy threshold, it has a small thermal 
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average at low temperatures. The magnitude of our thermal average cross section is of the 
same order as the meson exchange model calculation in ref. JT4]] with a cut-off A = 1 GeV. 
The shaded area in Fig. 7 give an evaluation of the uncertainties in our calculation due to 
the two procedures used to extract the hadronic amplitudes. 


V. SUMMARY AND CONCLUSIONS 


We have evaluated the hadronic amplitudes for the J/if dissociation by pions using the 
QCD sum rules based on a three-point function using vaccum-pion correlation functions. 
We have considered the OPE expansion up to twist-4 and we have worked in the soft-pion 
limit. Our work improves the former QCDSR calculation, done with a four-point function 
at the pion pole [^], since we have included more terms in the OPE expansion. We have 
shown that the twist-3 and twist-4 contributions to the sum rules are small when compared 
with the leading order contribution, showing a good “convergence” of the OPE expansion. 
We have checked that, taking the appropriated limit, we recover the previous result of [|22] . 

From the theoretical point of view, the use of QCDSR in this problem was responsible 
for real progress, being a step beyond models and beyond the previous leading twist calcu¬ 
lations § 0. i, n, m- This is specially true in the low energy region, close to the open 
charm production threshold. At higher energies our treatment is less reliable due to the 
approximations employed. 

Although a more sophisticated analysis of our uncertainties is still to be done, the shaded 
area in Fig. 6 shows that we can make some unambiguous statement concerning the behavior 
of <J n j/ip with the energy a/s. Our cross section grows monotonically with the c.m.s. energy 
but not as fast, near the thresholds, as it does in the calculations using meson exchange 
models [T3. [13], [L4], O, [16] . We have also shown the importance of respecting chiral symmetry, 


since the increase of the cross section near the threshold is strongly intensified when chiral 
symmetry is broken. In other words, our results suggest that, using meson exchange models 
is perfectly acceptable, provided that they include form factors and that they respect chiral 
symmetry. With these precautions, they can be a good tool to make predictions at somewhat 
higher energies. 

We have also evaluated the thermal average of the J/if — tt dissociation cross section. It 
increases with the temperature and at T = 150 MeV we get (a nJ ^v) ~ 0.2 — 0.4 mb which 
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is compatible with the values presented in Fig. 5 of ref. [|14|, i.e., in a meson exchange model 
with monopole form factor with cut-off A = 1 GeV. The use of this information will reduce 
the uncertainties in the calculations of the hadronic lifetimes of J/^, which are needed in 
simulations like those of ref. [Rfl. 
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APPENDIX A: SUM RULES FOR THE PROCESS ,J/'ip ix -> D D* 


Using r 3 = 7 ^ and in Eqs. (^|) and (j|), we obtain the following sum rules for 


the structures in Eq. flT5|): 

Af D * + A° d *M 2 + B\ )rr M 4 = 


Cdd* Idd*(M 2 ) 


+ /3tt 3 + 


2 M 2 


[-f*m c 

,-ml/M 2 

M 2 


(Al) 


A? B * + A? D *M 2 + BP d *M 4 = 


Cdd* f dd* (Yf 2 ) 


fnm c + 


4 (qq) 


m c (4:(qq)m c 5 2 f w ( 5 m 2 c ' 

'6M~A~17~ —{ 9 + Tp, 


3 /„ 

- m 2 /M 2 

M 2 ’ 


(A 2 ) 


\P D * + a° d *m 2 + b° d *m 4 = 


Cdd* f dd* (AT 2 


2 f n m, 


m r 


( 8{qq)m c S 2 f n ( ^ 10m 


6M 2 V fn 


M 2 


8 {qq) 
3 /tt 

—m 2 /M 2 


M 2 


(A3) 


Af D * + A° d *M 2 + B» d *M 4 = 


m 2 c (qq) e m C M2 

Cdd* f D D*(M 2 )Trr M 2 ’ 


(A4) 
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A? D * + A? D *M 2 + B° d *M a = 


( qq) e m2 / M2 


CdD* fDD*(M 2 ) fir M 2 


(A5) 


where 


foD*{M 2 ) — 


m ^ — m 2 D * 


e ~m 2 D /M 2 _ e ~m^/M 2 e ~m 2 D /M 2 _ e ~m 2 Dt /M 2 

~ m D m i>* — m D 


and 


Cdd* — 


m. 


m 2 D m D *m^f D f D * fy 


(A6) 


(A7) 


AP- D * and B? D * are the parameters introduced to account for double pole-continuum and 
single pole-continuum transitions respectively. 


APPENDIX B: SUM RULES FOR THE PROCESS J/i/j vr -*• D D 


Using r 3 = i'f 5 and T 4 = in Eqs. (||) and (|8|), we obtain the following sum rule for the 
structure in Eq. (|l(|) p4|: 


Add + AddM 2 + 5nnM 4 


m 2 e- m ' /M2 
M 2 


e -m 2 D /M 2 e -m 2 D /M 2 _ g -m^/M 2 


M 2 


- mf) 


j 2 m c (qq) fj 2 ^ | 5m 2 ' 


3 f n M 2 18 M 2 


M 2 


(Bl) 


APPENDIX C: SUM RULES FOR THE PROCESS J/i/j n ^ D* D* 


Using r 3 = 7l/ and T 4 = in Eqs. (||) and (||), we obtain the following sum rules for the 
structures in Eq. ©• 


Af • D ' + Af u 'M‘ + B"~"~ M' = 


D*D* 7\ , tjD* D* /I ,r4 


+ 


Cd*d* fD*D*{M 2 ) 

f n b 2 ( 25 5 m 2 \ 


-fir + 


2 m c (qq) 

3 UM 2 


6 M 2 \ 3 


+ 


3 M 2 I 


,-ml/M 2 

M 2 


(Cl) 


A 


+ A°* D * M 2 + B?* D * M* = 


Cd*d* /d*d*(A7 2 ) 


f 


f^b 2 /13 5m 2 

6M 2 l y + UP 


2 m c (qq) 
T 3 f n M 2 

—m 2 /M 2 


M 2 


(C2) 
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A D* D * + A D*D* m 2 + B D*D* m 4 = 


+ 


Cd*d* Id*d*(M 2 ) 


~ U + 


2m c (qq) 

SUM 2 


fJ 2 ( 

6 A/ 2 l 3 ' 3M 2 


19 5 m 2 ' 

+ 


D -m 2 /M 2 

M 2 


(C3) 


a?* d * + aU d *m 2 + b?* d *m 4 


Cd*d* fD*D*(M 2 ) 

fj 2 [23 5m 


2/, 


4m c (qq) 

SUM 2 


e ~m 2 c /M 2 


3M 2 \3 3M 2 / 


M 2 


(C4) 


Af D * + Af *■°*M 2 + B°* d *M a = 


Cd*d* f d*d* (Af 2 ) 


— /vr + 


23f n 5 2 

18 M 2 


,-ml/M 2 

M 2 


, (C5) 


Af - D * + Af D * A/ 2 + Bq* d * M 4 = 


Cd*d* /d*d*(A/ 2 ) 


~ /vr + 


25f n 5 2 
18 M 2 


o-ml/M 2 

M 2 


, (C6) 


A?* d * + AU d *M 2 + b?* d *m 4 


Cd*d* f d*d* (Af 2 ) 


/.- 


2m c (qq) 

SUM 2 


Sm 

M 2 


D -ml/M 2 

M 2 


(C7) 


a?‘ d * + aU d *m 2 + bF d *m 4 = 


Cd*d* f d*d*{M 2 ) 


u 


2 m c (qq) 
SUM 2 


2 mjto, 

M 2 


£—m 2 /M 2 

M 2 


(C8) 


A?‘ D * + A°* d *M 2 + bF d *m 4 


Cd*d * fD*D* (Af 2 ) 
13m 2 5 m 4 ^ 


^c-Ar , ^c(gg) 


M 2 / 53 

3 y 6 6M 2 ' 12M 4 / 


2 

e ~ m c/M 2 


+ 


2 U 


A/ 2 


(C9) 


A? 0 * D ‘ + A? 0 * D * M 2 + B® 0 * D * M 4 = 


Cd*d* f d*d*{M 2 ) 


2 m c (qq) 

3 U 


+ rn c f 37l 


£—m 2 /M 2 

M 4 


, (CIO) 


Aft 12 + A" " M 2 + 5" M 4 = 


m c (qq) e m c/ M2 


C D * D * f D *D*(M 2 ) U M 2 ’ 


(Cll) 


Af 2 * D * + Af 2 * D * A/ 2 + By2°* M 4 = 


Cd*d* f d*d* (AT 2 ) 


2m c / 371 


/X 


,-m 2 /M 2 

AT 4 


, (C12) 
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A?; d ’ + A°* D *M 2 + Dy^°* M a = 


f n 8 2 e m c/ M2 


Cd*d* f d*d*(M 2 ) 3 M A 


where 


fo*D*{M 2 ) — 


m D * — m 


i> . 


-m 2 D JM 2 e -m 2 D JM 2 _ p -m*/M‘ 


e d 


e w 


M 2 


~ m D* 


and 


C 


D*D* - —2 


™n*™*pfD*U 


(C13) 


(C14) 


(CIS) 


APPENDIX D: SUM RULES FOR THE MESON DECAY CONSTANTS 

For consistency we use in our analysis the QCDSR expressions for the decay constants of 
the J/A, D* and D mesons up to dimension four in lowest order of a s : 


_ 3 m c r*» (s — m c ) (ml _ s)/M 2 

}d 8tt 2 m 4 D U s 


rrc 


rri 


(qq)e 


( m 2 D —m 2 )/M 2 


(Dl) 


D 


fh = 


r s D* 


87 j 2 m 2 D * 


ds 


m 


2\2 


2 + ^ 
, s 


x g (m 2 D ,-s)/M 2 


m r 


m 


■(qq)e { 


(m 2 t -m 2 )/M 2 


D* 


^ 47T 2 Jim 2 


(s + 2m 2 ) yjs — 4 ?n 2 


^ ^ ^" Vcj y ' D ^'"‘ c c (m 2 -s)/M 2 


s 3/2 


(D2) 

(D3) 


where Sm stands for the continuum threshold of the meson M, which we parametrize as 
Sm = (rriM + A s ) 2 - The values of Sm are, in general, extracted from the two-point function 
sum rules for f D and /A* and /</, in Eqs. ([Dip , ( [I ) 2 | ) and m- Using the Borel region 
3 < < 6 GeV 2 for the D* and D mesons and 6 < M^ < 12 GeV 2 for the J/A, we found 

good stability for /A, /A* and /p with A s ~ 0.6 GeV. We obtained fo = 160 ± 5MeV, 
f D * = 220 ± lOMeV and = 280 ± lOMeV, which are compatible with the numerical 


values in Eq. (|T8|). 
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